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Aim To describe the organization, field work, forensic an- 
thropological examination, and DNA analysis conducted to 
identify the victims from a World War II mass grave found 
on the Dalmatian island of Daksa near Dubrovnik (Croatia) 
in 2009. 

Methods Excavation of the site was performed accord- 
ing to standard archeological procedures. Basic anthropo- 
logical examination was made to determine the minimum 
number of victims, sex, age at death, and height.The bones 
with pathological and traumatic changes were identified. 
DNA was extracted from powdered bones and relatives' 
blood samples. Y-chromosome and autosomal short tan- 
dem repeats (STR) were used to establish the relationship 
of the remains with the putative family members. 

Results The remains were found to belong to at least 53 
distinctive victims. All were male, mostly with gunshot 
wounds to the head. DNA analysis and cross-matching of 
the samples with relatives resulted in 1 4 positive identifica- 
tions using the Y-chromosomal STRs and 4 positive identi- 
fications using the autosomal STRs. 

ConclusionsThis study showed that even in cases of more 
than 50-year-old, highly degraded human remains from 
mass graves, Y-chromosomal and autosomal STRs analysis 
can contribute to identification of the victims. 
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After the partisan (communist-led) forces had entered Du- 
brovnik on October 18, 1944, mass executions followed, 
among them the execution on the nearby island of Daksa 
(Figure 1 and 2). According to historical sources (1 ), the ex- 
ecution took place on October 25, 1 944 and the following 
days. According to the verdict issued by the "Court of the 
Military Command for the South Dalmatian Region" and 
announced on a poster on October 29, 1 944, 36 Dubrovnik 
citizens, most of them prominent intellectuals, were "sen- 
tenced to be shot by a firing squad." The same historical 
sources do not state whether a trial was held or if there was 
a court under such a name (1). 



Figure 1 . 




Map of Croatia showing its position in Europe and Dubrovnik with its 
nearby island of Daksa. 

Figure 2. 




Map of the island of Daksa showing two sites of the mass grave, Loca- 
tion 1 and Location 2. 



In forensics, if the recovery of skeletons takes places a 
long time after death the possibilities to obtain informa- 
tion about victims'clothing, personal belongings, sex, age, 
height, and other characteristics may often be reduced. 
In such cases, DNA typing techniques may provide useful 
information. Y-short tandem repeat (STR) markers are pa- 
ternally inherited and allow identification of male missing 
persons when the only available reference is a male pater- 
nal relative. The autosomal STR DNA profiling was estab- 
lished by using the commercially available miniSTR assay, 
which proved useful to obtain information from ancient 
and degraded DNA samples and was used in cases when 
the only available reference was a missing person's daugh- 
ter (2). The DNA technology, including both STR analysis 
and mitochondrial DNA analysis, was already confirmed as 
a method of choice in the identification of missing persons 
in the 1 991 -1 995 war in Croatia (3,4). 

The aim of this study was describe the organization, field 
work, forensic anthropological examination, and DNA 
analysis by chromosome and autosomal STRs conducted 
to identify the victims from a World War II mass grave on 
the island of Daksa. 

MATERIALS AND METHODS 

The field work took place over fewer than 5 working days. 
During that time (including the time of departure/return 
from the island), the site was marked and trial trenches 
were dug and two locations treated (surface larger than 80 
m 2 , more than 1 50 m 3 of the excavated material) (Figure 3). 
Approximately 1 0 000 bones and skeletal fragments, and a 
number of items were discovered, packed, and transport- 
ed (Figure 4). In such circumstances, in most cases it was 
impossible to determine with certainty the body position 
and their interrelations. 

Field work 

The exhumation of mortal remains on the island of Daksa 
started on September 24, 2009, pursuant to the Order of 
the County Court in Dubrovnik. According to the historical 
data (1 ), it was assumed that the location of the mass grave 
was in the ruins of the former farm building basement. An- 
other mass grave was also believed to be in the immediate 
vicinity. The second location was discovered about 15 m 
away from the first location. 

Location 1 was an area surrounded by remains of a wall. It 
was divided, labeled, and staked out by sticks and ropes 
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into 55 equal fields measuring approximately 1 x 1 m. 
The surface layer of the soil was removed by a bulldozer. 
Bones started to emerge at the depth of approximately 
1 m from the surface (Figure 5). Due to a large number 
of mixed and separated bones in different positions, in- 
termingled with large stones, it was decided to give up 
the skeleton presentation in situ. This primarily referred to 
the Location 1 , fields labeled "d.'The bones were carefully 
excavated. After having defined the orientation (cranial/ 
caudal) and the position related to the marked fields and 
the layer depth, the bones were packed in labeled card- 
board boxes. 

Location 2 was not staked out as Location 1 , but was shaped 
as the excavation went on. Three parallel trenches were 
dug and marked as 2, 3, and 4, and divided into 9 equal 
fields measuring approximately 1 x 1 m. The trenches fol- 
lowed the direction of cracks in limestone. The findings of 
the excavation were the same as those at Location I.The 
work was additionally obstructed by thick pine roots. Vari- 
ous items were found along with the bones at both loca- 
tions: fragments of priestly collars, rosaries, buttons, neck- 



Figure 3. 




Mass grave site. 
Figure 4. 




laces, bullets, and shells (Figure 6). The excavation went on 
up to the hard layer at the depth of about 2 m. 

Work at the pathology room 

The remains were transferred to the Department of Pathol- 
ogy of the Dubrovnik General Hospital, where they were 
cleaned and sorted out. Where possible, the bones were 
placed in anatomical position and analyzed to determine 
the age, sex, height, and changes that had occurred during 
lifetime or after death, and photographed. Samples for the 
DNA analysis were taken subsequently. After having taken 
measures of thighbones and head diameters, the height of 
a person was determined by theTrotter equation for calcu- 
lating the body height (5). 

Figure 5. 




Human skeletons in the mass grave after the surface layer of the soil was 
removed by a bulldozer. 

Figure 6. 




Mass grave exhumation. 



The objects found at the site. 
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TABLE 1 . Archeological data from the mass grave at Daksa Island, Location 1 . The data are divided according to marked fields, see 
Figure 7 

Field Scapula Clavicula Humerus Ulna Radius Ribs Sternum Vertebrae Femur Pelvis 1/2 Sacrum Tibia Fibula Small bones Skull 



Da 


2 


2 


2 


2 


2 


20 


1 


20 


1 


1 


0 


3 


4 


82 


1 


Ea 


3 


3 


3 


4 


3 


29 


2 


36 


5 


5 


3 


3 


3 


155 


1 


Fa 


3 


3 


3 


2 


3 


33 


1 


32 


0 


1 


1 


0 


0 


80 


2 


Ga 


0 


0 


0 


0 


0 


1 


0 


4 


2 


0 


0 


4 


4 


97 


2 


Ha 


5 


6 


6 


6 


6 


84 


3 


59 


6 


6 


3 


2 


2 


176 


0 


la 


0 


0 


0 


0 


0 


0 


0 


6 


2 


0 


0 


3 


4 


68 


1 


Ja 


2 


2 


2 


1 


1 


19 


1 


28 


1 


4 


2 


0 


0 


63 


1 


Ka 


2 


2 


2 


2 


2 


17 


1 


19 


0 


0 


0 


0 


0 


26 


1 


lb 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


2 


] 


15 


o 


Jb 


0 


0 


0 


1 


1 


0 


0 


0 


1 


0 


0 


1 


0 


0 


0 


Total 


17 


18 


18 


18 


18 


203 


9 


204 


18 


17 


9 


18 


18 


762 


9 


Fc 


0 


0 


0 


0 


0 


7 


0 


7 


4 


3 


1 


5 


6 


0 


0 


Gc 


4 


4 


3 


1 


1 


31 


1 


23 


1 


3 


2 


1 


1 


81 


1 


He 


7 


4 


6 


6 


6 


52 


5 


63 


10 


6 


6 


3 


3 


200 


1 


Ic 


3 


3 


5 


3 


4 


33 


2 


33 


1 


4 


2 


0 


0 


93 


2 


Kc 


5 


4 


5 


6 


5 


75 


2 


58 


0 


3 


3 


0 


0 


188 


2 


Fd 


0 


0 


0 


0 


0 


0 


0 


1 


4 


1 


1 


3 


3 


13 


0 


Gd 


4 


4 


6 


7 


6 


81 


4 


58 


8 


7 


1 


14 


14 


281 


3 


Hd 


4 


5 


3 


3 


5 


59 


2 


69 


3 


4 


2 


4 


4 


175 


3 


Id 


7 


7 


7 


9 


8 


60 


2 


71 


3 


3 


1 


12 


13 


428 


5 


Jd 


4 


4 


3 


2 


3 


68 


2 


51 


12 


6 


4 


4 


2 


86 


2 


Kd 


6 


7 


8 


6 


8 


62 


3 


65 


0 


3 


0 


0 


0 


173 


4 


Total 


44 


42 


46 


43 


46 


528 


23 


499 


46 


43 


23 


46 


46 


1718 


23 


Overall 


61 


60 


64 


61 


64 


731 


32 


703 


64 


60 


32 


64 


64 


2480 


32 



TABLE 2. Archeological data from the mass grave at Daksa Island, Location 2. The data are divided according to marked fields (No. 
2-5), see Figure 8 

Field Scapula Clavicula Humerus Ulna Radius Ribs Sternum Vertebrae Femur Pelvis 1/2 Sacrum Tibia Fibula Small bones Skull 



2 


2 


3 


4 


1 


3 


23 


1 


31 


3 


3 


2 


2 


2 


25 


2 


2a 


2 


3 


3 


4 


4 


35 


2 


43 


5 


5 


2 


5 


4 


75 


0 


2b 


2 


2 


1 


2 


2 


27 


1 


23 


0 


0 


0 


4 


6 


231 


2 


2c 


2 


0 


2 


1 


1 


4 


0 


13 


4 


10 


5 


7 


6 


28 


0 


2d 


3 


3 


3 


3 


3 


57 


2 


55 


6 


1 


0 


0 


0 


35 


2 


2e 


6 


5 


3 


3 


3 


37 


3 


35 


0 


0 


1 


0 


0 


76 


2 


2f 


1 


3 


3 


4 


4 


27 


1 


34 


2 


1 


0 


2 


2 


89 


2 


Total 


18 


19 


19 


18 


20 


210 


10 


234 


20 


20 


10 


20 


20 


559 


10 


3 


2 


3 


2 


2 


2 


14 


0 


21 


2 


2 


1 


6 


8 


181 


3 


3a 


4 


3 


3 


4 


3 


37 


2 


40 


7 


2 


1 


6 


4 


69 


0 


3b 


2 


1 


3 


3 


4 


36 


2 


28 


5 


10 


4 


1 


1 


115 


1 


3c 


6 


7 


6 


3 


3 


68 


3 


66 


0 


0 


1 


1 


0 


57 


2 


3d 


0 


0 


0 


2 


2 


4 


0 


9 


0 


0 


0 


0 


1 


17 


1 


Total 


14 


14 


14 


14 


14 


159 


7 


164 


14 


14 


7 


14 


14 


439 


7 


4 


6 


6 


5 


3 


4 


69 


3 


36 


4 


5 


2 


2 


1 


67 


2 


4a 


0 


0 


0 


1 


2 


0 


0 


6 


2 


0 


0 


4 


5 


84 


1 


Total 


6 


6 


5 


4 


6 


69 


3 


42 


6 


5 


2 


6 


6 


151 


3 


5 


0 


0 


1 


1 


0 


0 


0 


0 


0 


0 


0 


1 


0 


3 


1 


Overall 


38 


39 


39 


37 


40 


438 


20 


440 


40 


39 


19 


41 


40 


1152 


21 
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Boundary values for the sex determination, based on the 
long bone length, were determined as described by Slaus 
et al (6). Where possible, sex was determined by pelvic 
bones and skull characteristics. Due to the significant level 
of damage, this was generally not possible. Also, because 
of a large number of bones in a limited area and in the 
same layer, it was impossible to determine with certainty 
the body parts that belonged to the same person. 



in which they were found (eg, Location 1 field Da). Fields 
in which skeletal remains were not found are not listed in 
the tables. 

We made a schematic presentation of the assumed posi- 
tion of a particular body in relation to other bodies and 
their spatial orientation within the labeled schematic site 
(Figure 7 and 8). 



The age was determined by the pelvic articular surface 
(symphysis) following the Suchey Brooks method (7). Since 
most of the bones were badly damaged, a reliable deter- 
mination was generally not possible. The location and the 
area of discovery and the total number of the bones dis- 
covered are shown in Table 1 and Table 2. Skeletal remains 
were identified and classified according to macro-sites (Lo- 
cation 1 and 2) and micro-sites (marked fields-quadrants) 

Figure 7. 




DNA analysis 

Prior to the DNA extraction, the bone/tooth surface was 
cleaned by abrasion with a grinding tip and sandpaper, the 
bone/tooth was crushed into small fragments, and stored 
in sterile polypropylene tubes at -20°C until analysis. After 
this initial treatment, the DNA extraction was performed 
as described by Alonso et al (8). At least two independent 
extractions were done (one or two craniums that present- 
ed teeth were sampled, as well as different postcranial 
bones like femurs). The DNA from blood and bloodstain 
reference samples of the living relatives was extracted by 
standard Chelex 100 protocols (9). The polymerase chain 
reaction (PCR) amplifications were performed on Parkin- 
Elmer thermal Cycler 9700, using the AmpFiSTR Yfiler PCR 
Amplification Kit (Applied Biosystems, Foster City, CA, USA) 
and MiniFiler PCR Amplification Kit (Applied Biosystems) 
(1 0,1 1 ). Typing of PCR products was performed on the ABI 



Figure 





Schematic presentation of the finds at Location 1 . 



Schematic presentation of the finds at Location 2. 
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Prism 310 Genetic Analyzer (Applied Biosystems) with the 
Data Collection Software. Electropherogram data were an- 
alyzed with the GeneMapper ID Software, version 3.2 (Ap- 
plied Biosystems). The internal standard was Liz-500 (Ap- 
plied Biosystems). 

Analysis of typing results 

DNA profiles from skeletal remains were analyzed and 
compared with DNA profiles of the living relatives. The 
database was kept in the Microsoft Access 2000 (Micro- 
soft, Seattle, WA, USA). Microsoft Excel 2000 (Microsoft) 
was used for statistical calculation. Calculation for statisti- 
cal probability of biological relationship was performed 



Figure 9. 




Gunshot wounds in the occipital region of a victim. 

Figure 10. 



GENERATION 




Schematic presentation of positive identification in a case where great- 
grandfather's brother was identified by matching with great-grandson's 
paternally inherited Y chromosome. 



according to standard protocols and the data from the Y 
Chromosome Haplotype Reference Database (YHRD) da- 
tabase (2,12,13). 

RESULTS 

Human remains were found at two sites designated as 
Location 1 and Location 2, and the skeletal remains be- 
longed to at least 53 persons (Table 1 and 2). It can be stat- 
ed with certainty that the bones were human and more 
than 50 years old, the time corresponding to the World 
War II. The skeletal remains had distinctive male anthropo- 
morphic characteristics. Due to the bone damage of the 
remains, the age of a certain number of persons could not 
be determined. 

Many of them had metal fillings, gold, and silver teeth and 
dental bridges and dentures. 

Gunshot wounds were identified in 22 persons. Practi- 
cally all identified entry gunshot wounds were localized 
in the occipital region (Figure 9). One gunshot wound in 
the occipital region was identified in 1 5 persons, 2 gunshot 
wounds in 5 persons, and 3 gunshot wounds in 2 persons. 
Only one gunshot wound was identified in an area other 
than the skull, ie, the hip socket. Many skulls were found in 
fragments with peri-mortem fractures, although the head 
gunshot wound was not found. A tibia bone with a peri 
mortem fracture was also found on the site. 

DNA was successfully extracted from 49 out of 53 skeletal 
remains (Table 3). In each case of tooth and bone samples, 
several DNA extractions were performed. AmpFiSTR Yfiler 
PCR Amplification Kit produced 34 out of 49 profiles, 32 of 
which were full profiles and 2 partial profiles. MiniFiler PCR 
Amplification Kit produced 40 out of 49 DNA profiles and 
only 2 were partial. Identifier PCR Amplification Kit was 
used for 4 skeletons, and 2 full and 2 partial profiles were 

TABLE 3. Summary of DNA analysis data for the skeletal re- 



mains from the island of Daksa 
Characteristics 

Exhumed remains 53 

Obtained DNA profiles 49 

Ratio of obtained profiles (%) 92.43 

Relatives' samples 23 

Positive identifications 18 

Ratio of positive identifications (%) 33.96 

Ratio of positive identifications (%) regarding 78.26 
the number of requesting families 
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Y chromosome short tandem repeat loci mutations in a father/son pair 
with double peaks at DYS437 locus and with null allele at DYS448 locus 
(in circles). 

obtained. DNA was also extracted from 23 blood samples 
of the living relatives. 

Cross-matching of the results with the relatives' data re- 
sulted in 18 positive identifications. In 4 cases, the victims 
were identified by using only the autosomal STR systems 
for identification because positively identified remains be- 
longed to the fathers of the daughters who gave the blood 
sample. The results obtained by MiniFiler PCR Amplification 
Kit were confirmed by the Identifier PCR Amplification Kit, 
which allows simultaneous amplification of 15 autosomal 
STR loci, whereas MiniFiler PCR Amplification Kit gives the 
profiling result for 8 autosomal loci (both kits include sex- 
determining locus amelogenin). 

In other 14 cases, the victims were identified by the Y- 
chromosomal STR system because the victims were pa- 
ternally related to the individuals who gave blood sam- 
ples. Nine positive identifications were made for father/ 
son pairs, 3 for nephew/uncle pairs, 1 for grandchild/ 
grandfather pair, and 1 for great-grandson/great-grand- 
father's brother pair (Figure 10). We also found rare Y 
chromosome STR loci mutations in a father/son pair with 
double peaks at DYS437 locus, as well as with the null al- 



lele at DYS448 locus observed in the same father/son pair 
(Figure 1 1). 

DISCUSSION 

This study presents the anthropological, forensic, and ge- 
netic analysis of victims from the mass grave on the island 
of Daksa, near Dubrovnik. The remains were found to be- 
long to at least 53 distinctive victims. All were male, most- 
ly with gunshot wounds to the head. DNA analysis and 
cross-matching of the samples with relatives resulted in 
14 positive identifications using the Y-chromosomal STRs 
and 4 positive identifications using the autosomal STRs. Fo- 
rensic analysis revealed that most persons with identified 
gunshot wounds to the head died a violent death due to 
those wounds. The same applies to the persons with iden- 
tified peri-mortem skull fractures and fragmented skulls. 
Arm bones and fore-arm bones tied up with wire indicate 
that this was a mass execution grave. The manner of the 
burial and particularly the localization of the penetrating 
gunshot wounds of 9 mm in diameter and the 9 mm pis- 
tol bullets with corresponding shells found at the site also 
suggest that this was an execution. The persons executed 
were probably in the kneeling position at the burial site or 
in its immediate vicinity. The objects found at the site indi- 
cate that they were mostly civilians. Also, according to the 
items found in the mass grave and to the results of DNA 
identification, at least 3 victims were priests. 

One of the main problems in excavations of the mass 
graves dating from the World War II is a large number of 
damaged and separated skeletal remains (14). The exca- 
vations are usually complicated by environmental factors, 
such as difficult access to the site, subsequently raised con- 
structions, vegetation, and often only assumed micro-lo- 
cation of the gravesite (15). Financing, manpower issues, 
the available time, and particularly the lack of experienced 
and educated teams also have a significant impact on the 
methodology of excavation. 

The main objective of the World War II mass graves exhu- 
mation is to determine the minimum number of victims (if 
not always the precise), identify the cause and manner of 
death, and carry out the identification test, if possible, and 
return the remains to their families, in accordance with the 
Geneva conventions (16). 

Since the bones of a large number of victims were found 
separated and intermingled, often in piles, a fast and 
bulk approach is satisfactory, particularly if the grid 
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use (network) and the orientation of long bones at the site 
is carefully observed. Such method provides a significant 
saving of time and resources and still gives necessary an- 
swers to the questions raised. With the methods described 
above and additional work in the pathology room and at 
the desk, it is possible to subsequently reconstruct posi- 
tions and relations of the remains by combining the data 
collected during the exhumation with the basic anthropo- 
metrics data revealed in the pathology room. 

Anthropometrics data and the teeth status have an ex- 
tremely limited role in the identification of mass graves vic- 
tims dating from the World War II, mostly because there are 
no ante-mortem data related to victims (1 7,1 8). There was a 
relatively limited number of citizens of Dubrovnik who were 
expected to be found on the island of Daksa, but after more 
than 60 years no relevant data on their physical characteris- 
tics (more accurate than short/tall) or the teeth status could 
have been obtained, although a large number of victims 
had significant prosthetics works. There were no significant 
problems in the anthropological determination of the vic- 
tims' sex, whereas the age determination (especially for the 
older ones) was not sufficiently precise for identification. We 
were able to determine the actual age of the victims only af- 
ter the DNA analysis and positive identification, according to 
the information given by the relatives (19). 

It was easy to identify the cause and the manner of death, 
since the execution with a single shot to the back of the 
head was a common manner of execution in mass graves 
victims from the World War II discovered in this region, as 
was the burial at the place of the execution (1 7). 

A well-organized field work and careful proceedings in the 
pathologist room can "tell the story"aboutthe last moments 
of the victims' lives, while the identification of victims is exclu- 
sively in the domain of the DNA analysis (14,1 7,20). The lat- 
ter was the imperative in our research, and was mostly done 
through Y chromosome analysis of male skeletal remains. 

We successfully identified 1 4 persons by the AmpFiSTR Yfil- 
er PCR Amplification Kit. One of them was an honorable 
Croatian priest, reverend Petar Perica, known as the author 
of some of the most beautiful Croatian religious songs (Box 
1 ). His identification was confirmed by the Y-chromosome 
STR profile compatibility with the blood sample of his 
brother's great-grandson. 

While observing a null Y-STR locus should not repre- 
sent a problem for the profile interpretation, more 



BOX 1 . Croatian religious poem written by Rev. Petar Perica, 
one of the identified victims 

Heavenly Virgin, Queen of the Croats 
(Rajska djevo kraljice Hrvata) 

Hail, Mary, full of all graces, 
Eternal sunshine, clad in radiance, 
Circling your brow, a starry crown, 
Below your feet, hell's dragon groans. 
Heavenly Virgin, Queen of the Croats, 
Our Mother, Our Golden Dawn, 
From devoted hearts, receive a gift, 
Receive our pure and fervent love. 
Blessed are you, all Immaculate, 
The serpent's breath taints not your breast! 
Star of happiness, resplendent for us, 
Disperse evil darkness, nights of sin! 

http://www.youtube.com/watch?v=l_SC8_SGHvO 



than one peak at one or more loci could lead to a grave 
misinterpretation of a profile and its sample source. The 
same rare mutations and their interpretation have been al- 
ready reported (21). Multiple variations (deletion or dupli- 
cations), even if rare, can happen in the same single-source 
sample and these examples should be a reminder before 
drawing premature conclusions. Forensic interpretation of 
Y-haplotype profiles should be in the focus because mul- 
tiple alleles at various loci do not necessarily indicate that 
the sample originates from a mixture. This also showed the 
importance of the DNA identification methods and their 
application in the case work. 

This study showed that even in cases of more than 50-year- 
old, highly degraded human remains from mass graves, Y- 
chromosomal and autosomal STRs analysis can contribute 
to the identification of the victims. 
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